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We use Cosmic Microwave Background (CMB) observations together with the Red-sequence 
Cluster Survey (RCS) weak lensing results to derive constraints on a range of cosmological pa- 
rameters. This particular choice of observations is motivated by their robust physical interpreta- 
tion and complementarity. Our combined analysis, including a weak nucleosynthesis constraint, 
yields accurate determinations of a number of parameters including the amplitude of fluctuations 
as = 0.89 ± 0.05 and matter density Q m = 0.30 ± 0.03. We also find a value for the Hubble pa- 
rameter of Ho = 70 ± 3 km s _1 Mpc _1 , in good agreement with the Hubble Space Telescope (HST) 
key-project result. We conclude that the combination of CMB and weak lensing data provides some 
of the most powerful constraints available in cosmology today. 
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The physics behind the anisotropics we see in the mi- 
crowave background is well studied and understood 0. 
The evolution of the photon distribution function in the 
tight coupling era and through decoupling is well inside 
the linear perturbation regime and is the reason for the 
CMB's unique status as a probe of cosmological models. 
The physical interpretation of the angular power spec- 
trum of primary CMB anisotropics is unambiguous when 
restricted to the inflationary paradigm and given a suit- 
ably parametrized spectrum of initial perturbations. 

The recently released WMAP first year results Q have 
revealed the CMB angular power spectrum with unprece- 
dented accuracy to multipoles below I = 900 Q. The 
results are a stunning confirmation of the acoustic oscil- 
lation picture, with perturbations arising from an initial 
super-horizon spectrum of predominantly adiabatic fluc- 
tuations, as predicted for example by simple inflationary 
models. The measurements of the first two acoustic peaks 
has confirmed in precise detail earlier detections of the 
peak/dip pattern on scales below the sound horizon at 
last scattering H 0, 0, II H E3, EH 

On its own, the current picture of the CMB made up 
of the WMAP results together with high resolution CBI 
and ACBAR 1 observations implies tight constraints on 
a number of parameters |l2|, ll2j ; the curvature in units 
of critical density VLk and various other parameters in 
the combinations determined by the physical mechanisms 
which give rise to the observed CMB anisotropy. In ad- 
dition the measurement of a cross-correlation between 
the polarization and temperature anisotropy EH 0] is 
the first significant detection of reionization in the CMB, 
which gives a constraint on the optical depth to the last 
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scattering surface. 

Although the CMB data alone provide tight con- 
straints on some parameter combinations, other combi- 
nations are very poorly constrained due to partial degen- 
eracies. The addition of other data such as measurements 
of the matter power spectrum P(k) is essential to break 
these degeneracies and tightly constrain the parameters 
of most interest individually. One way to infer the matter 
power spectrum is to rely on visible tracers of the (dark) 
matter distribution such as galaxy redshift surveys or ob- 
servations of the Lyman-a forest. The Lyman-a forest 
gives a way to measure the linear power spectrum of neu- 
tral gas at redshifts higher than those probed by galaxy 
surveys. 

The combination of CMB, 2dFGRS [II], and Lyman- 
a forest data 

El 

yields tight constraints on the density 
of dark matter and vacuum energy, and also reveal an 
indication of a running of the scalar spectral index char- 
acterized by the parameter dn s /d\nk [T3 | . However in- 
ferring the matter power spectrum using these techniques 
involves a heuristic treatment of the relation between the 
tracers and the dark matter usually referred to as 'bias- 
ing' E3] • As we enter the much heralded era of precision 
observations, such heuristic treatments might limit the 
accuracy with which parameters can be determined. A 
direct measurement of the power spectrum would not 
suffer from such limitations. 

In terms of physical interpretation, measurements of 
the lensing signal induced by the LSS (cosmic shear) hold 
a unique position in the growing set of observational tools 
available to cosmologists; it is a direct probe of the pro- 
jected matter power spectrum over a redshift range deter- 
mined by the lensed sources and over scales ranging from 
the linear to non-linear regime. The intervening LSS in- 
duces a small, coherent correlation in the shapes of the 
background galaxies which nowadays can be measured 
accurately 

Ea E3, El S E3 The use of weak lensing 

data is not without challenges: the small signal requires 
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large survey areas and a careful removal of the observa- 
tional distortions. However separation of the shear sig- 
nal into gradient ("E-Type") and curl ("B-Type") com- 
ponents provides a control on systematics including the 
presence of intrinsic alignments of nearby galaxies or sys- 
tematically induced distortions in the image. The RCS 
53 sq. deg. results used in this work [24| have a low 
B-Type component on large scales together with a well 
determined redshift distribution of background sources. 

In this letter we present results from cosmological pa- 
rameter fits using only CMB and weak lensing data. The 
motivation for this approach is to provide constraints on 
parameters using only observables with robust physical 
interpretations. 

To evaluate the posterior distribution of the param- 
eters of interest from the data we use an extension of 
the publicly available Markov Chain Monte Carlo pack- 
age cosmomc 2 , as described in j^- We calculate the 
likelihood of each cosmological model with respect to a 
combination of CMB and RCS data. The CMB data con- 
sists of WMAP data below i = 900 and CBI , ACBAR , 
and VSA band powers above i = 800 where the WMAP 
data is noise dominated and hence the band powers are 
essentially independent. To compare each angular power 
spectrym to the WMAP data we use the likelihood calcu- 
lation routine made available by the WMAP team 3 [liij ]. 

For each model we also calculate the mass aperture 
variance ( M 2 (9 )) [2fij | at each aperture 9 sampled by the 
RCS results |24j . The mass aperture variance is a narrow 
filter of the convergence power spectrum P K (£) defined as 



P K (t) = 




(1) 

where \ is the radial coordinate and /kIx) * s the comov- 
ing angular diameter distance to x- P3D (k, x(-z))is the 3D 
power spectrum of matter fluctuations. For each model 
we use the matter power spectrum calculated by CAMB 
|29] | at z — and rescale to z > using the solution for 
growth of linear perturbations. To include the non-linear 
contribution to the power spectrum at each redshift we 
use the HALOFIT procedure |27j. The procedure has been 
calibrated using numerical simulations of structure for- 
mation and is significantly more accurate than the previ- 
ous procedure by Peacock & Dodds [2^] . In particular it 
reproduces accurately, with rms errors of a few percent, 
the full non-linear spectrum in standard ACDM mod- 
els down to scales k ~ 10/iMpc -1 . The accuracy of the 
HALOFIT procedure is adequate for current weak lensing 
data although future surveys will require more accurate 
estimates of the full, non-linear power spectrum. This 
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will most probably require the use of large numbers of 
numerical simulations to calibrate directly the non-linear 
evolution in the full parameter space. 

The function g( X ) = J* H d X 'p(x')fK(x' ~ x)/Sk(x') 
is the source-averaged distance ratio where p(x{z)) de- 
scribes the redshift distribution of sources in the shear 
survey which is approximated by the function p(z) ~ 
(z/z s ) a exp [-(z/z s f]. The values a = 4.7, /3 = 1.7, 
and z s — 0.302 give the best fit to the observed redshift 
distribution. To allow for the uncertainty in the mean 
redshift of the distribution we marginalize over the range 
of values z s £ [0.274,0.337] for each likelihood evalua- 
tion. This corresponds to the ±3<r range indicated by 
the x 2 °f the fit to the photometric redshift distribu- 
tion. The mean redshift for this choice of parameters is 
(z) = 0.54 — 0.66. We assume a Gaussian prior for z s in 
this range. 

For each model sampled by the Monte Carlo chain we 
calculate the log likelihood with respect to the RCS data 

hxL = -\ (<M|). - (M 2 ap )) Cr> ((i5f \ . - (M 2 ap)] ) , 

(2) 

where (M% p ). is the observed mass aperture variance at 
an aperture 9i and Cjj is the covariance matrix of the 
data [24| . This result is added to the log likelihoods from 
the CMB fit for the same model to obtain the full likeli- 
hood with respect to both CMB and RCS data. 

We sample the probabilities with respect to six ba- 
sic cosmological parameters: the physical densities of 
baryons flbh 2 , and cold dark matter £l c h 2 , the Hubble pa- 
rameter Hq = lQOh km s _1 Mpc _1 , a reionization redshift 
parameter z ro , and a constant spectral index n s and am- 
plitude A s of the initial scalar curvature perturbations. 
We assume the universe is spatially flat, with purely adi- 
abatic perturbations evolving according to General Rela- 
tivity. The density of a cosmological constant type com- 
ponent fl\ follows from £!a = 1.0 — f2 m . We generated 
sixteen converged Monte Carlo chains using the CMB 
data only, removed burn in and thinned to obtain fairly 
independent samples. The matter power spectrum and 
RCS likelihood was then computed for each sample, and 
importance sampling used to adjust the chain weights ac- 
cordingly (see 25] ) . The resultant set of weighted sam- 
ples for the full posterior distribution from the CMB 
and RCS data were then used to compute our results. 
The only external prior assumed is a conservative big 
bang nucleosynthesis (BBN) Gaussian prior of the form 
n b h 2 = 0.022 ± 0.002 (la) We include this prior to 
partially break the remaining n s -fli,h 2 -T~A s degeneracy 
in the CMB data. The action of this is similar to the 
t < 0.3 prior adopted in the WMAP analysis [HEH. 

From the set of samples it is simple to also compute the 
posterior distribution of other derivable quantities such 
as the rms amplitude of matter fluctuations on 8/i _1 Mpc 
scales assuming linear evolution, erg, the total matter den- 
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FIG. 1: One dimensional, marginalized probability distribu- 
tions for a selection of parameters. The dashed (black) curves 
are for CMB only. The solid (red) curves include the RCS 
data. We see the weak lensing data, while being consistent 
with the CMB only results, narrow down a number of distri- 
butions considerably. In particular the combination of data 
provides a tight independent constraint on Hq. The matter 
density and fluctuation amplitude are also much better con- 
strained with the combination of CMB and weak lensing than 
with just CMB. 



sity, f2 m , the optical depth to last scattering, r, and the 
age of the universe. In this letter we do not consider ten- 
sor perturbations, dynamical dark energy candidates, or 
a running spectral index. We will explore the constraints 
on these generalized models from CMB and weak lensing 
data future work. 

The set of samples from the full six dimensional pa- 
rameter space can be used to evaluate marginalized pa- 
rameter distributions by evaluating the weighted number 
density of samples ignoring the values of the parameters 
marginalized over. In Fig. 2] we show the one dimen- 
sional marginalized distributions for a number of param- 
eters. Each panel compares the distribution obtained us- 
ing CMB data with that obtained using CMB and RCS 
data together; both also include the weak BBN prior dis- 
cussed above. The effect of adding the weak lensing re- 
sults is clearly seen in a number of parameters. 

In Table [I] we summarize the marginalized constraints 
for a number of fundamental and derived parameters. 
We show the results obtained with and without inclusion 
of the RCS data. The addition of RCS data reduces 
the errors on cr s , Q m , H , Q\, and Q c h 2 . We also show 
constraints on the 'classical' combinations probed by LSS 
data, namely, the constrained direction <j$£l^ and the 
shape parameter T w Q, m h. 



TABLE I: Marginalized constraints for a selection of param- 
eters. The left column uses only CMB data, the right column 
is for CMB and RCS data. The only external prior included 
for both results is a Gaussian BBN prior of 0.022 ± 0.002. All 
errors are 68% confidence intervals. 
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It is instructive to look at the marginalized, two di- 
mensional likelihood in the {Sl m , cs) plane to understand 
how drastic improvements in the determination of the 
two parameters are obtained (Fig. 0). The RCS data 
alone is near degenerate in a particular direction while 
CMB data alone provides broad constraints in a quasi- 
orthogonal direction to RCS. The combination of the two 
data sets give a much tighter confidence region. The re- 
gion of intersection in six dimensions has slightly above 
average CMB likelihood, as is readily assessed using the 
importance weighted samples, so the data sets are highly 
consistent even in the full parameter space. The spread 
of the CMB posterior in the direction of the RCS de- 
generacy is largely due to the uncertainty remaining in 
the optical depth, as is clear for the tight constraint for 
a%e~ T given in Table [I] 

Overall our results are consistent with similar con- 
straints from a combination of CMB, 2dFGRS, and 
Lyman-a data [illlij] with similar or smaller errors. The 
values obtained for eg using the WMAP data are some- 
what higher than those obtained previously from CMB 
data due to the new evidence for a significant optical 
depth and a slightly higher anisotropy amplitude than 
previous observations indicated Q. This is still lower, 
although not inconsistent, with estimates of as from a 
possible Sunyaev-Zeldovich Effect (SZE) contribution to 
the CMB power spectrum at high-£. The latest estimates 
using the CBI deep-field results [UHi and ACBAR and 
BIMA .33] data suggests a value of a§ z = 0.98±£^ H3 
with large errors due mainly to the non-Gaussian nature 
of the SZE. Increasingly accurate measurements of the 
CMB power spectrum at high-^ will reduce these errors 
drastically and comparing the two independent determi- 
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much promise for future investigations. 
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FIG. 2: The two dimensional, marginalized likelihoods for 
the (fi m ,£78) plane. The overlaid, filled contours show the 
68% and 95% integration levels for the distributions. Bottom 
- RCS only, Middle - CMB only, Top - CMB+RCS. 



nations of a% will be useful in increasing our understand- 
ing of the cluster properties that determine the SZE. 

Our result for the Hubble parameter is consistent 
within la with the HST key-project result [3^ but has 
smaller errors. Similarly, the value for the matter density 
fl m = 0.30±0.03 is consistent with other determinations. 
The addition of RCS data leave estimates of the scalar 
spectral tilt n s essentially unaffected. This is due to the 
small range of scales probed by the RCS weak lensing 
results. Future surveys will most certainly have much 
more leverage on n s as they will probe a range in scales 
an order of magnitude larger. 

We have shown how CMB and weak lensing results 
can be combined to obtain constraints on cosmological 
parameters that depend on observations that have sim- 
ple physical interpretations. Although only first gener- 
ation weak lensing data are currently available our ap- 
proach yields results with errors comparable to or even 
smaller than those obtained using CMB in combination 
with other types of surveys. These results are encourag- 
ing for the use of next generation weak lensing surveys in 
deriving robust parameter fits. In particular the Canada- 
France-Hawaii- Telescope (CFHT) Legacy Survey ~ 170 
sq. deg. cosmic shear project will be a major step for- 
ward in the field of weak lensing. 

The increasing accuracy in the determination of the 
source redshift distribution in future surveys will also 
help in reducing uncertainties and systematics tied to any 
intrinsic alignment in the ellipticity of nearby sources. It 
will also introduce the possibility of resolving separate 
redshift contributions to the convergence power spectrum 
(Eq.G) thus enhancing the parameter fitting ability of the 
observations. 

We conclude that the combined CMB, weak lensing ap- 
proach to parameter determination already constitutes a 
competitive alternative to other combinations and holds 
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